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It has been proposed that a quantum computer can be constructed based on electron spins 
in quantum dots or based on a superconducting nanocircuit. During two-qubit operations, 
the fluctuation of the coupling parameters is a critical factor. One source of such fluctuation 
is the stirring of the background charges. We focused on the influence of this fluctuation 
on a coupled spin qubit system. The induced fluctuation in exchange coupling changes the 
amount of entanglement, fidelity, and purity. In our previous study, the background charge 
fluctuations were found to be an important channel of dephasing for a single Josephson qubit. 
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Among the various proposals for quantum computation, quantum bits (qubits) in solid 
state materials, such as superconducting Josephson junctions 1 and quantum dots, 2-4 have the 
advantage of scalability. Proposals to implement a quantum computer using superconduct- 
ing nanocircuits are proving to be very promising, 5-9 and several experiments have already 
highlighted the quantum properties of these devices. 10-12 Such a coherent- two-level system 
constitutes a qubit and the quantum computation can be carried out as the unitary opera- 
tion functioning on the multiple qubit system. Essentially, this quantum coherence must be 
maintained during computation. However, dephasing is hard to avoid due to the system's 
interaction with the environment. The dephasing is characterized by the dephasing time T2, 
and various environments can cause dephasing. 

Background charge fluctuations (BCFs) have been observed in various kinds of sys- 
tems. 13-16 In nanoscale systems, BCFs are electrostatic potential fluctuations arising due 
to the dynamics of an electron, or hole, on a charge trap. In particular, the charges at charge 
traps fluctuate with the Lorentzian spectrum form, which is called random telegraph noise 
in the time domain. 16 ' 17 The random distribution of the positions of such dynamical charge 
traps and their time constants leads to BCFs or l//noise. 18 In solid-state charge qubits, these 
BCFs result in a dynamical electrostatic disturbance and hence, dephasing. It should be noted 
that this dephasing process does not mean the qubit being entangled with the environment, 

*E-mail address: toshifumiitakura@scphys.kyoto-u.ac.jp 



1/10 



J. Phys. Soc. Jpn. Letter 

but rather, that the stochastically evolution of an external classical field is suppressing the 
density matrix elements of the qubit after averaging out over statistically distributed samples. 

We had shown that BCFs are important channel of dephasing for a single Josephson charge 
qubit system. 19 ' 20 In the present study, we investigate the effect of BCFs on the two-qubit- 
gate operation. To construct a controllable quantum computer, one requires the suppression of 
dephasing and accurate universal quantum gate which consists of single qubit operations and 
two-qubit operations. Therefore, to address these manipulations, we examine the dephasing 
of coupled qubit system, which is experimentally current topic and is urgent to analyze what 
cause of dephasing is important in these systems. 

Recently, it has been shown that the interaction between electron spin in a quantum 
dot and environments is weak, 21 ' 22 and one can expect a very long dephasing time of an 
electron spin. Therefore, the proposal of quantum computer using electron spin in a quantum 
dot is promising. 4 For the electron spin qubit, however, the effects of ; the fluctuation of local 
magnetic field, 4,23 ' 24 and the spin-orbit interaction 24 are important. Moreover, the fluctuation 
in the exchange coupling is important during the gate operation, because, for two-qubit gate 
operation, one uses the exchange interaction between the two quantum dots. This type of 
dephasing occurs when the charge of traps change the coulomb energy of two-qubit state and 
fluctuate the exchange energy between the two qubits. Therefore, we investigate the effect of 
BCFs on the two-qubit gate operation. 

We examine the time evolution of two-qubit density matrix which obeys the two-qubit 
Hamiltonian H = J(i)Si ■ S2, where J(t) depends on time, and we neglect the magnetic field 
on each qubit system. We examine the amount of fluctuation in exchange coupling. (Figl.) 

We define the quantum variables of two qubits' coordinate as f{ and r~2, and define 
the environment variables of a charge trap and its near-by electron reservoir as 7*3 and 
7*4. We define the Hamiltonian of two qubit system as double well in a two-dimensional 
layer (z=0) such that H = Ei=i,2 ^ + C,h = ^p 2 + V{n),C = ^L^,V(x,y) = 
( 2Z y-^5-(x 2 — a 2 ) 2 + m ^-y 2 ), where m is the effective mass of an electron in the quan- 
tum dot, lo the confinement frequency of confinement potential, and 2a the distance be- 
tween two potential minima. The exchange energy of electron spins between the dots is 
given by J = (c{e~ d2 h{d 2 )) + f (1 + d 2 )) where c = ^/2(e 2 /a B )/hw, I (x) is 0- 

th order Bessel function, and d = a/ 'as, clb = y^h/muj. 25 The interaction Hamiltonian is 
V\ = e 1 1 \J\f\ — i%\ 2 — e 2 / \J\f2 — fl\ 2 ^ e 2 /r + e 2 ((r~i +7*2) ■ (r~3 — fX))/r z , where r ~ y^r"^ 2 ~ 
y^rll 2 . From the above calculation, the dynamic part of the interaction Hamiltonian be- 
tween charge trap and qubit system is given by V^fi,^) = e 2 (^+^)^3-^)+fai+i/2)fe-i/4) ? 
where we set z\ and zi to zero. Using the Heitler-London approximation, we define, \S) 
and \T) to be singlet and triplet states such that, \S) = (|12) + |21))/ v / 2(l + S 2 ), |T) = 
(|12) - |21))/ v / 2(l - S 2 ), where, S 2 = |(1|2)| 2 = exp(-2d 2 ) is an overlap integral and 
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Fig. 1. Scheme of coupled qubit system and charge 
trap 
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aligns in y-direction, the fluctuation of the exchange coupling does not exist. Then, the amount 
of fluctuation due to a single charge trap is written by, J lsing \ e = ^$3?) llK ^^fe^ ) 2 - 25 
For many traps, the magnitude of fluctuation in the exchange coupling becomes as follows. 
The dipole vector of BCFs is ^(sin^jcos^sintfjsin^cos^) at distance r^, where p is dis- 
tance between a trap and electron reservoir. Then the magnitude of fluctuation is given by 

r _ V J, . , _ V 13 1 / e 2 psin^cos^q ^2 1 2ir e 4 1 Nja 2 p 2 w u„ r „ at. :„ t u„ 

J\ ~ Z^lsmglc - l^i^-dpy n 5 > sinh(2^) ~ ~3~ foU sinh(2d 2 ) ~r^dF~ ' Wllere iV * 1S tile 

density of charge traps and is the distance where the dipole approximation becomes invalid 
(a < r d ). 

Next, we examine the time evolution of qubit density matrix. The time evolution of the 
qubit system is defined by the unitary operator 

l7(t)= e -*/oJMSi-s a *-_ (!) 



The exchange coupling operator is expressed by using permutation operator P\i such that 

P12 + I 



Si • S2 



P\2 has following properties: 
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where m is integer. We assume the fluctuation in the BCF obeys the random telegraph-type 
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Fig. 2. Scheme of Bloch sphere with 
pointing moving quantum state. 

(0, 0),a=(0, 0),/?=(0, f0,7=(7r/2, 0),6=(n/2, ir),e={n/2, tt/2). 

noise. J(t) = Jo + J\X(t), where X(t) takes the value of 1 or with characteristic time r. 
The density matrix at time t is given by p{t) = (U(t)p(t = 0)U*(t)}, where U(t) is unitary 
operator which describes the time evolution in terms of the bases | jj), | ij), | ||), | jj) and 
() is the ensemble average about the stochastic process. 20 ' 26 ' 2 To accomplish two-qubit gate 
operation, the restricted subspace spanned by | jj), I IT) i s enough because other states are 
decoherence free from this type of dephasing.In restricted sub-Hilbert space, the initial pure 
state is given by |^) = cos|| ||) + sin|e _l?i | ||), and the quantum state is represented by 
using a Bloch sphere as shown in Fig. 2. We obtain the density matrix at time t after taking 
the ensemble average and calculate T2 as had been done for a single qubit system. 20,26 ' 27 The 
states on one-dimensional line that connects two maximally entangled states (|(| ||) + | J.|)) 
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and II) — | II)). do not evolve during the gate operation. The distance from the origin 
of the Bloch sphere to quantum state represents the purity. The time evolution trace of 
the quantum state changes from sphere to ellipsoidal because of dephasing. The maximally 
entangled states |^) = (| ||) =F I do not evolve in time, since they are eigenstates 

of the Hamiltonian. Fig. 3 (a),(b) shows time dependence of quantum state with some choice 
of initial condition. We estimate T 2 _1 for many charge traps with different characteristic time 

as follows. We assume the transition times between occupied and empty states are equal and 

-w 

the temperature dependence of transition time obeys the thermal activation type r = Ae kBT , 
where W is thermal activation energy and A is characteristic time scale which is independent 
of temperature. To estimate the effect of many charge traps, we average over the magnitude of 
fluctuation in the exchange coupling and thermal activation energy. Since the dephasing rate 
by a single trap is J 2 t7j for weak coupling and H/2t for strong coupling, 19 ' 20 total T 2 _1 is given 
by T 2 ~ 1 = 3 2W) ^ l where 1/Wo is distribution of thermal activation energy. We define the gate 
operation time as t p = tt / Jq. Then, the gate quality factor is given by Q = 2~2/t p ~ 3 ^ . 

Next, we examine the quantum information quantities of the coupled qubit system. First, 
we examine the criterion of entanglement. This criterion comes from the negativity of mini- 
mum eigenvalue of a partially transposed qubit density matrix. 28 ' 29 Fig. 4 shows time depen- 
dence of criterion of entanglement. For Fig. 4-6 we set T2/T = 2. Fig. 5 shows time dependence 
of fidelity. The minimum eigenvalue oscillates with time. When sin 9 cos <fi = 0, separable states, 
ID) I II)) appear between entangled states during time evolution. Whereas, the qubits are 
always entangled if sin 9 cos <fr 7^ 0. Next, we study the fidelity of p(t) relative to p(t = 0), 
defined by J~(t) = Tr[p(t)p(0)], which shows accuracy of the quantum gate. 30 ' 31 Starting from 
the maximally entangled state, fidelity is 1. For more general initial conditions, fidelity is given 
by 

T = + cos (i J t)e~ t/T2 cos 2 9 



h 



+ 



- cos(i J t)e~'/ T2 sin 2 9 sin 2 <f> 

_■ lb 



+ i sin 2 9 cos 2 (f). (4) 

Finally, we examine the purity (V = Tr(p(t) 2 )), which is related to linear entropy of qubit 
system as Su n = 1 — p. 30 ' 32 The analytical expression of purity is 



1 1 2 
V = — I — sin 9 cos 

2 2 



2 



+ i(sin 2 #sin 2 (/> + cos 2 #)e- 2 */ T2 . (5) 

Figure 6 shows the time dependence of purity. The purity decreases monotonically, which 
means that the entropy of the qubit increases due to the dephasing. The initially pure qubit 
state becomes a mixed state, namely, V(t = 0) = 1 becomes V(t = 00) = \ + \ sin 2 0cos 2 <p. 
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Fig. 3. Time depnedence of quantum state, (a) initial 
state a=(9 = 0,(f) = 0). (b) initial state 13— (9 = 
0/4,^ = tt/2). 
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Fig. 4. Criterion of entanglement. 
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Fig. 6. Time dependence of purity. 
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